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Abstract 

We discuss constraints on which flat directions can have large vacuum expectation values 
(VEVs) after inflation. We show that only flat directions which are not charged under B-L 
and develop positive pressure due to renormalization group efi^ects can have large VEVs 
of order Mp\. For example, within the MSSM only the H^Hd flat direction is found to 
be viable. This strongly constrains the embedding of a broad class of inflationary models 
in the MSSM or some other supersymmetric extension of the SM. For flat directions with 
negative pressure, the condensate fragments into very large Q-balls which we call Q-giants. 
We discuss the formation, evolution and reheating of these Q-giants and show that they 
decay too late. The analysis requires taking into account new phases of the flat directions, 
which have been overlooked in the formation and dynamics of the Q-balls. These constraints 
may be ameliorated by invoking a short period of thermal inflation. The latter, however, is 
viable in a very narrow window of parameter space and requires flne tuning. 



1 Introduction 



Inflation is a successful cosmological paradigm both from an observational point of view and 
theoretical point of view. It is observationally successful as the mounting cosmological data 
indicate that we live in a FRW universe with a flat spatial section and that density perturbations 
that seeded structure formation is adiabatic and nearly scale invariant. It is even partially 
falsifiable as a large class of inflationary models called "slow-roll" inflationary models may 
be ruled out through the consistency condition PP if we are successful in measuring B-mode 
polarization which contains information about primordial gravity waves. From a theoretical 
point of view, inflation is appealing because it is the only known way to honestly solve the 
flatness problem of cosmology. Furthermore, because the energy scale of its physics may be 
as high as 10^^ GeV, fundamental theories such as string theory may have a good chance of 
making contact with its observable physics. 

Unfortunately, inflation does not have any smoking gun signatures yet because it has 
not made any connection with controlled experimental settings such as terrestrial colliders. 
Furthermore, there are no compellingly beautiful models of inflation as the constraints from 
the requisite large number of e-foldings, the small amplitude of density perturbations, and the 
smallness of the departure of the spectral index from scale invariance necessitate unnaturally 
tuned parameters in ad hoc constructions of inflationary models. The best hope for discovering 
whether inflation has occurred in our universe's history seems to be in flnding a natural inflaton 
candidate in physics beyond the standard model (SM) discovered at future colliders. 

One of the best motivated ideas for physics beyond the SM is supersymmetry (SUSY), and 
in particular, signatures of the minimal supersymmetric standard model (MSSM) are actively 
sought after in the next generation of high energy physics experiments such as the LHC (for 
a recent review in this direction, see for example, j^). Furthermore, engineering required to 
construct successful inflationary scenarios is often more natural in the context of supersymme- 
try. Among the plethora of scalars afforded by supersymmetrizing the SM, flat directions are 
particularly enticing candidates for inflatons.^ Finally, from a perturbative string theory point 
of view, one of the simpler embeddings of the SM into the MSSM arises from gravity mediated 
SUSY breaking context. Hence, we are naturally led to ask which flat directions in the gravity 
mediated SUSY breaking scenario of MSSM may be part of the inflaton dynamics. 

To this end, the story of Q-balls, a class of nontopological solitons [BlISISl) plays a crucial 
^Throughout this paper, we collectively denote flat directions by 0. 
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role. Because nearly all of the flat directions in the MSSM have an accidental U{1) symmetry at 
the renormalizable operator level, they fulfill one of the conditions for the existence of Q-balls 
(that of having a conserved global charge) [HI E] . The other crucial condition for the existence 
of Q-ball is that the scalar field potential U satisfies 



mm 



U 



< lu'\o). (1) 



This condition also has a good chance of being satisfied for MSSM flat directions because its 
interactions with the heavy degrees of freedom renormalize the mass (which is of the order of 
the SUSY breaking gravitino mass ^T^3/2) of the flat direction as 1 + Kln{jp-) and if 

K <0, Eq. (HI) can be satisfied (Sj.^ 

This story finds relevance in cosmology [HI because for situations in which a flat direction 
is displaced far from the origin at the end of inflation, one can show that when K < 0, there is 
negative pressure llOj, which generates the necessary instability for Q-ball formation. Further- 
more, since the U{1) is an accidental symmetry, non-renormalizable operators which break this 
U (1) induce large charge in the condensate which ends up inside the Q-balls. 

Therefore, as is known |51l^ lll[ll2U13[m] . Q-balls generically form at the end of inflation 
if a flat direction is turned on during inflation and in particular, if it plays a role in the 
inflationary dynamics. It is also known that if the Q-balls are too large and survive certain 
violent destructive processes, they will live a long time before decaying, leading to a reheating 
temperature too low to be compatible with cosmology. Hence, certain MSSM flat directions 
can be ruled out as candidates for the inflaton (or more precisely, flat directions cannot be 
displaced far from the origin at the end of inflation) if the Q-balls form along these directions 
and live a long time. 

We show that Q-balls with VEVs of order the Planck scale (which we call "Q-giants" due to 

their attendant large charge) generically live a long time, thus eliminating many flat directions 

from the possibility of being displaced far from the minimum at the end of inflation. The new 

physical observation that allows us to make this conclusion is that between the Q-giants, the 

density of gas of flat direction particles is necessarily large, thereby inducing large masses to all 

particles such that the decay rate of the gas is highly suppressed. Consequently, not only do the 

Q-giants live longer due to their suppressed decay rates, but also the thermal plasma produced 

by the decay of the gas is too cold to dissociate the Q-giants. Thus, all MSSM flat directions 

with negative value of K are excluded as inflaton candidates. The only viable flat direction of 

^Integrating out many heavy degrees of freedom during inflation is generic since the MSSM flat directions can 
be displaced far from the minimum, inducing large masses to the flelds to which it couple. 
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the MSSM that does not carry B-L charge is found to be H^Hd- Finally, although B-L charge 
density carried by the effective Affleck-Dine field with VEV close to Mp\ is generically too 
large to be consistent with cosmology even when no Q-giants form, such directions may survive 
due to charge dilution through mechanisms such as thermal inflation (although as we discuss, 
constructing such thermal inflationary models will be contrived and requires flne tuning). Thus, 
HuL direction which has positive pressure may survive as an inflaton candidate. 

Although most of the paper is concerned with the MSSM flat directions having Mpi fleld 
variations (hence, from a string theory point of view, these would naturally be moduli fields), 
it is important to note that most of the results can be applied to other systems in which there 
is Q-giant formation starting from a field condensate. In particular, the constraints hold for 
flat directions in any extension of the SM, such as GUT or string embedding. Indeed, there are 
realizations of inflation where some of the moduli fields cannot be part of the MSSM, and should 
therefore arise from flat directions of some extension ^21- The results of this paper may strongly 
constrain a given extension. Furthermore, during (F-term) inflation, flat directions generically 
obtain corrections which may drive them alway from the origin |161 117j . The existence of such 
directions, even without playing a role in inflation, is highly constrained in the same manner. 

The outline of the paper is the following. Section [21 contains an executive summary of 
the cosmological scenario and the results. Section |31 contains a discussion of the framework 
and some of the models of inflation to which it applies. Sections and [S] discuss the case of 
negative pressure fiat directions. Section |1] deals with the formation of Q-giants after infiation 
and section[21deals with their subsequent evolution leading to too low of a reheat temperature to 
be compatible with cosmology. In sectionlHlwe discuss positive pressure flat directions and show 
that they are cosmologically acceptable. Section [7| deals with thermal inflation as a possible 
loophole to our results. We argue that such a scenario however requires fine tuning and is very 
restricted. We conclude in section |H1 

We will use the notation Mpi = l/^/SvrG/v = 2.4 x 10^^ GeV throughout the paper. 

2 The Cosmological Scenario: An Executive Summary 

In this section, we briefiy outline the cosmological scenario and the main results of this paper. 

We study infiationary models in the framework of gravity mediated supersymmetry break- 
ing. In many natural models of infiation an inflationary epoch is induced with the use of moduli 
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whose VEVs are Mp\ away from the true minimum and which dominate the energy density dur- 
ing inflation. A canonical example is hybrid inflation |18 ( I15 ( ITO l 120 1 [2 H I22j . Such moduli are 
usually related to flat directions in a moduli space of a supersymmetric theory. The purpose of 
this paper is to investigate the dynamics of such flat directions and thus place constraints on 
their possible role in inflation. Furthermore, in many theories flat directions are driven away 
from the origin due to Hubble- induced couplings |161 117j . The existence of such flat directions 
may again be constrained. We concentrate on the MSSM, but stress that the analysis holds for 
any supersymmetric model with gravity mediated SUSY breaking. 

2.1 Assumptions 

Let us spell out the basic assumptions for which our scenario holds: 

• At the end of inflation, the energy density is dominated by a flat direction, </>, which is 
completely flat in the supersymmetric limit, 7713/2 — > 0. As a consequence its potential is 
proportional to w-o/a and is therefore of the form. 



• One consequence of the above potential is that (j) generically begins its evolution at Mpy. 

• Another consequence is an approximate U{1) for {(j)) <C Mpi. Thus at the onset of its 
evolution (j) is "side-kicked" such that it carries an order one charge which is later released 
to the light particles including those of the MSSM. An immediate conclusion is that (j) 
cannot carry B-L charge jl9j with reheating temperature above the sphaleron transition 
temperature except under certain exceptional situations.^ 

2.2 Cosmological Scenario 

Let us now point out the main features of the cosmological evolution that follows from the 
above assumptions. After taking into account renormalization group (RG) corrections to the 
modulus potential, the latter takes the form [51 lllj. 




(2) 



U{(t>)=ml 1 + Klog 



M2 




) 



k=l 



00 




+ C.C. 



(3) 



^More on this point and its caveats in section [7| 
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where ~ ""^3/2) ^ is the scale at which 772-3/2 is defined, typically the GUT or Planck scale 
and K is related to the renormalization group beta function. K can be negative or positive and 
typically ranges between 10~^ to 10~^. The values of depend on the symmetries (possibly 
discrete R-symmetries) of the model. Initially the condensate oscillates and due to the RG 
corrections positive or negative pressure is built up, depending on the sign of K 10 . Before 
we present our computation, we summarize the (j) field evolution in each of the two cases as 
follows: 

Negative Pressure, K <{) 

• Negative pressure fragments the coherent condensate into very large Q-balls (Q-giants) 
and non-coherent, non-thermal particles with typical momentum of order rn^^ . 

• While almost all of the initial charge is stored in Q-giants, both the Q-giants and the (j) 
gas initially store comparable amount of energy. 

• Q-giant collisions break the giants into either smaller Q-balls or into the non-coherent 
gas. We show that while the gas comes to dominate the energy density rather quickly, at 
least a percent is left in the Q-giants. 

• Outside the Q-giants (</)) = 0, but 

(^'f^) ~ lQ~'^\K\^Ml^. (4) 

This finite density effect induces large masses to any fields which are directly coupled to 
(f). As cj) decays through theses heavy particles, its decay rate is suppressed by l/{4>^4>). 

• The gas reheats the light particles outside the Q-giants to intermediate temperature, 

Tintcr ^ 10' GeV. (5) 

This temperature is too low to dissociate the Q-giants, and the smaller than naively 
expected temperature is due to the decay rate being suppressed by 

• At this stage the temperature outside (for example, of the MSSM) dominates the energy 
density. However, since the Q-giants still hold at least a percent of the energy density in 
the form of non-relativistic matter, there is sufficient time for them to overcome the ther- 
mal radiation and dominate the energy density a long time before they decay, reheating 
to, 

Treheat ^ 10"^ GcV. (6) 
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Positive Pressure 

• In such cases, no fragmentation occurs. 

• Due to the positive pressure, a finite Jeans scale develops, preventing the growth of small 
perturbations. 

• By the time the evolution of the perturbations become non-linear, the Jeans scale is large 
enough to prevent the formation of primordial black holes (BHs). Therefore, there is no 
significant cosmological constraint coming from overproduction of BHs. 

2.3 Conclusions 

The above dynamics strongly constrains the possible use of flat directions in inflation: 

• Negative pressure fiat directions cannot develop VEVs of order Mpi for viable 
cosmology. They are excluded since they decay very late, dominating the energy density 
during nucleosynthesis thus destroying its predictions. Furthermore, if charged under B, 
they release too much baryon number which cannot be erased at such low temperature. 

• Positive pressure flat directions do not have the above problems, and further solve the 
problem of extensive formation of primordial BHs. 

• Without invoking thermal inflation [23L 1241 125| . B-L directions that either carry B or 
which reheat to temperatures high enough to produce B from sphaleron transitions are 
excluded due to excess baryon production. In section [7| we stress the difficulties in writing 
a model which incorporates a period of thermal inflation. 

• For the case of flat directions in the MSSM, the only positive pressure flat 
direction which is not charged under B-L is HuH^- Hence within the MSSM and 
excluding exceptional situations presented in section [3 this is the only possible direction 
which can play a role in hybrid inflation or any other modular inflation. 

• More generally, flat directions obtain Hubble-induced contributions to their masses during 
inflation. Such contributions drive them towards or away from the origin depending on 
the sign of the mass. Hence HuH^ is the only exactly flat direction in the MSSM which 
can obtain negative contributions of this kind. 
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3 The Framework 



One way of classifying models of inflation is by considering whether or not the inflaton field 
varies a field distance of order Mpi during inflation. This classification has phenomenological 
utility as it can tell us about whether significant gravitational waves are generated or not 
j26j . Theoretically, large field variations signal theories with symmetries as potentials must be 
extremely flat to generically allow such large field variations without costing excessive energy. 
For example, in the context of perturbative string theory, Planckian field variations are natural, 
as the most ubiquitous scale in this theory is the string scale, and supersymmetry and modular 
invariance protect flat potential directions. As we will see, Mpi field variations are also natural 
in hybrid inflationary models embedded in supergravity. 

The results of our present paper regarding the formation and longevity of Q-giants are 
likely to apply to most scenarios of inflation in which the inflaton field at the end of inflation 
must still travel Planckian field values to reach the current minimum, and if the potential 
for this field satisfies the condition for Q-ball formation. The reasons is of course because 
the defining characteristic of a Q-giant, large Q-charge, is due to the large field displacement 
from the potential minimum. Therefore, even though we devote most of our attention in this 
paper with the end of inflation in the context of MSSM flat directions having Planckian field 
displacements, in the next subsection, we discuss more generally possible inflationary models 
where our results might be applicable. Note that although non-renormalizable operators are 
generically important for the inflationary scenarios with (near-)Planck field variation, the Q- 
giant analysis is insensitive to the details of the non-renormalizable operators. 

3.1 Applicable models 

There are many models of inflation that have been proposed which can involve Planck scale 
field variations. Arguably, most popular models of such type are chaotic inflationary scenarios 
[27\ and hybrid inflationary scenarios I28j . Typically, for chaotic inflationary scenarios, the 
field variation is of order Mp\ during inflation and after inflation. On the other hand, hybrid 
inflation can involve field variations that are very small during inflation while the field variation 
is large after the end of inflation. 

In the context of SUSY, the generic structure of SUSY inflation comes from the scalar 
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potential 

U = lma:b']DaD, + e^/^^p.5*^*(Al^)(I?ii^)* - -J^e^/^p.VF*H^ (7) 
^ PI 

where 

9ij = d^,d^*K (8) 
dW 1 dK 

and where the D-terms in the flat space hmit look like 

Da--gia-g{<^*T^j<^j) (10) 

with fab being the gauge kinetic function and not necessarily vanishing for a U{\). For [/ 7^ 0, 
if Da 7^ 0, we have D-term inflation [291 12U1 OT] . while otherwise, we have F-term inflation for 
which Ty / or D^W / 0. 

F-term models typically suffer from the r] = U" MpjU > 1 problem (the 77 problem) |2H] 
which arises from the e^^^Pi factor. Although the D-term models typically do not suffer from 
this problem (since the F-term is zero by construction), in the context of string theory, if 
g£,a comes from an anomalous U{1), its value is generically too large to be compatible with 
cosmology. Note that even in D-term models the rj problem can arise since 

l^[f~^^]DaD,r^eg\l + ^ + ...) (11) 

which contributes 

V^- T^— (12) 

' FA2 Sir ^ ^ 

which would generically be too large for inflation to take place [30]. However, this problem can 

be cured by imposing appropriate discrete symmetries. 

D-term inflation typically realizes hybrid inflationary scenarios while F-term inflationary 
scenarios can realize both hybrid and chaotic scenarios. Because chaotic inflationary scenar- 
ios generically require the mass of the inflaton to be of order 10^^ GeV, these directions are 
generically heavier than MSSM flat directions (i.e. this scenario is largely irrelevant for embed- 
ding MSSM flat directions into an inflationary scenario). Furthermore, although this scenario 
naturally will have Q-ball production, it gives rise to sufficiently large reheating temperature 
for a viable cosmology. Hence, we will now focus on discussing the salient features of hybrid 
inflationary scenarios. 
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A common hybrid inflationary potential parameterization is 

= .^(m=-^)%^ + !I!^ + 2,. (13) 

where a is the field direction that is rolling during inflation when (p = and (p is the waterfall 
field direction that shuts off inflation at the end when a rolls down to o"^ < = 2M^/t^/A^. 
In the standard scenario, we have 

K^M"^ > mV^ > ja* (14) 

which implies that the constant energy term k?M^ dominates the energy density during in- 
flation. The dominance of the constant energy density is the defining characteristic of hybrid 
inflationary models. When combined with the slow roll conditions and number of e-fold con- 
straint with m = 771-3/2 ITeV and M = f3M-p\ (the latter says that field variations of the 
waterfall field (j) are Planckian at the end of inflation), the phenomenologically allowed param- 
eter ranges are 

- <^m^f« (15) 



"^3/2 



^K^P^ < 10-^9 (16) 



< ^2^4 ^ -{SttYPti ~ 10-« (17) 



^1/2 N 



,2^4 "'3/2 

ln(^ 

where P-ji ~ 10^^ is the adiabatic scalar perturbation power spectrum, ~ 50 is the number 
of e-folds required for inflation, and 

1 1 (k2m4)3/2 

aN ~ (19) 

is the field value N e-folds before the end of inflation. ^ The conditions given by Eqs. (|15() 
and H16() can be associated with Eq. ()14|) . while Eq. H17|) comes from the standard slow roll 
conditions e = ^{U'/Uf <^ 1 and rj = M^^{U" /U) < 1. 

A few remarks are now in order. 



1. In the absence of symmetries, the parameters k, A, and 7 must be finely tuned. 

*This is under the approximation that inflation really ends when a = cr/ which is only approximately true 
when the field is not too light. 
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2. Within the context of supergravity, the waterfall field, 0, is Mp\ away from its minimum 
at the end of inflation |T3]. This has two consequences: (i) (p is weakly coupled to the 
MSSM and therefore suffers from a moduli problem jl9j . One solution to this problem is 
to let the true minimum (as well as the origin of field space) be an enhanced symmetry 
point in which case the decay width of (p is large, (ii) By the time (j) decays, a has mass 
of order the Planck scale. 

3. If the parameterization of Eq. (|13() is used, the waterfall field (p rolls away from the origin 
at the end of inflation. As noted by 15_, this is one of the main reasons why it is difficult 
to identify one of the MSSM flat directions with 0: MSSM flat directions have VEVs 
that are zero or much smaller than Mpi today. One may naively think it is easy to shift 
the origin of the inflationary potential. However, this is incorrect since a vacuum shift 
generically breaks the underlying symmetries governing the form of the original potential. 
Nevertheless, it could well be, if one is to solve the moduli problem, that the light MSSM 
fields do live away from the origin on another enhanced symmetry point. While common in 
string theory, such theories are complicated to write down, since the symmetries cannot be 
realized linearly. One should therefore think about the above Lagrangian as an expansion 
of the full theory around the origin, and can thus view as an MSSM field rolling away 
from the origin, and into the MSSM enhanced symmetry point. 

As a concrete example in which is rolling away from the origin, let us first consider one 
of the the models of jj^ (for improvements on this type of model and further discussion see 
^ni)- They use the superpotential 

(where g' is a coupling constant) and an instanton induced soft breaking terms with scales which 
are natural for gravity mediation SUSY breaking scenario. The resulting potential (neglecting 
the SUGRA corrections) can be written approximately as 

2 2 2 

U{^, « - + + 1^ (I'^l VP + kh^n (21) 

which can be identified with Eq. (|13j) through k ~ am^/2/{Mpi)^ 7 ~ S'^l'/'P/^pii ^^'^ ^ ~ 
(7^|cjAr|^/-/Vfp[. Note that this choice of k with a ~ 1//3 is natural for a fiat direction field (p in 
supersymmetry. With this parameterization, because a is rolling towards smaller values during 
inflation, the condition ~ 10^/3^m3/2 ^ cj/ ~ 1^^3/2/^ must be satisfied. Hence, /3 cannot 
be too small to inflate sufficient number of e-folds, which means, as we advertised, that Planck 
scale field variation of (j) is natural in the context of SUSY hybrid inflation |19| . 
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For an example of a model in which the waterfall field ends up at i?i> = 0, consider the model 
of [Snj. This model has the merit in that it also solves the rj problem. The supersymmetric 
part of the model is defined by the superpotential 

W = AiaxVi + ^2j—^(p^i^2 (22) 

and a D-term 

D = A^- \x\^ - \cl)f + iViP + nlVsp. (23) 

If IxP + |0P ^ A'^) potential is minimized for ijji = '4>2 = 0. Hence, ignoring for the moment 
the soft terms, the inflationary scalar potential becomes 

U = A^kPlxl^ + A^SSl'/'r + ^(A^ - Ixl^ - \<t^\'? (24) 
which inflates while a is large enough 

H^''/^^(^)"^'- 

with X = 3,nd the waterfall field pinned at 

i*i^-^^-"|(]i|^)""*^- (^«> 

Here we have expanded to leading order in A2- 

During inflation, and taking into account soft breaking terms, the effective potential is of 
the form 

A 2n-4 

U = \l-^,K' + ml,a\ (27) 

In the parameterization of Eq. ()13() . this model can be identified with M ~ A and k ~ 
A2(A/Mpi)"'~^. Hence, to satisfy CMB fluctuation amplitude, one find using Eqs. (|T9|l and 
((23), A ~ Mpi and A2 — m3/2/Mpi, giving U ~ 7713^2^?! during inflation. 

After the end of inflation, we have x rolling out to A while (j) rolls to 0. Note however, that 
due to the D-term, cf) is not a flat direction of the full gauge group. Thus, if one is to embed the 
MSSM in this model, one must introduce an additional U{\) gauge group. As a final remark, 
we point out that the solution of j29j to the T/-problem involves having W = = which 
means that (j) is F-flat. Thus these class of models may incorporate flat directions of the MSSM 
or any other supersymmetric extension. 

In all of these models, there exists a global U{1) which is important for the formation of 
Q-balls as we later explain. We will see in section [3.21 these accidental U{1) symmetries are 
generic in MSSM flat directions. 
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3.2 A Specific Example - Flat Directions in the MSSM 



Flat directions are ubiquitous in supersymmetric theories. In particular there are numerous 
such directions in the MSSM, as was systematically found in |3H I17j . As discussed in the 
introduction, our goal is to constrain the embedding of moduli in models of inflation in the 
MSSM and more generally the existence of flat directions in models which involve an inflationary 
epoch. One such restriction was explained in ^M,- Indeed, if we are in a situation where all 
the energy in the universe after inflation is stored in some field with large VEV, then this 
field cannot carry B-L number, or else too much baryon asymmetry is generated. This kind of 
reasoning assumes, of course, no fine tuning,^ although the amount of this fine tuning is not 
clear once fiuctuations are taken into account, due to processes similar to tachyonic preheating. 
Thus we are left with only five possible fiat directions for the waterfall field: 





B + L 


HuHd 





QuLe 





QuQd 





QQQL 


+2 


uude 


-2 



Further constraints on the embedding may be pursued by investigating the cosmological 
evolution of such fiat directions. Indeed it is known that under certain conditions (i.e. negative 
pressure and such fiat directions being charged under some global U{\)), they fragment to form 
Q-balls. To see that such a U{\) exists for the above directions, we note the following. In 
the MSSM, the above directions are flat at the level of the renormalizable superpotential, and 
one expects these flat directions to be lifted by non-renormalizable terms. However, in the 
models discussed in section ITTl in order for the flat direction to be naturally M-p\ away from its 
minimum, it must be an exact fiat direction in the supersymmetric limit ^H]- This in particular 
means that soft breaking terms are proportional to rn^i^ |^. Hence assuming for simplicity 
only one fiat direction, denoted by (/>, its potential is therefore, 

^(0) = ml/^Ml, U{r/M^„ |</.|'). (28) 

^Finely tuned initial conditions could easily allow for an eccentric orbit which would suppress the produced 
baryon number. 
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Thus as long as (j) does not couple to light fields,^ it has an approximate U{1) for (^/Afpi < 1. 
For the directions not charged under B or L, this is no longer true once other fields become 
light. However, as we will see in section 14.31 with the exception of HuH^, this only occurs 
very late in the cosmological evolution, due to (i?5>^</>) condensate. For the H^H^ direction the 
above does not hold, due to the ^-term which gives a mass to (j) already in the supersymmetric 
limit. Here we do not attempt to solve the /x-problem and assume either is of order m-3/2 or 
it is promoted to a field as in the NMSSM [^. Due to the special properties of this direction 
(positive pressure), the absence of an approximate U{1) is irrelevant and as we show below it 
is most suitable to play the roll of the modulus in inflation. 

A remark is now in order. Above we have discussed the approximate U{1) in the scalar 
sector. One might worry that outside Q-balls, (p may decay into fermions which remain light. 
However, as we show below by appealing to supersymmetry, since bosons which interact directly 
with (j) become heavy, the fermions also becomes heavy and therefore the Q-balls remain stable. 
Hence below we assume that along the flat direction the interaction to lighter SM fields (to 
which it can decay) is through some particles whose mass is proportional to the VEV of (j) or 



Flat directions are constrained even without driving inflation. Within the context of su- 
pergravity, flat directions obtain Hubble-induced corrections to their masses, 



during inflation |161 117j . Here Hj is the Hubble constant at inflation and c is a model-dependent 
dimensionless parameter. If c > 0, (/> becomes tachyonic during inflation and is quickly driven 
away from zero to a minimum which, for exact flat directions, is Mpi away. It then remains at 
this minimum until the end of inflation, when it starts rolling back towards the origin. This 
situation is constrained by the analysis below, much like the embedding of modular inflation. 

4 Negative Pressure Flat Directions - Formation of Q-Giants 

This section is devoted to analyzing the matter content right after the condensate fragments. As 

we show, the universe contains both Q-giants and a non-thermal distribution of low momenta 

(j) particle.^ This will serve as an initial configuration for the system which we evolve through 

the reheating period in section |3 

^Light fields means that these fields are not integrated out. 
^This is in contrast with the thermal effects discussed in |33II34| . 




(29) 
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In the case of gravity mediated SUSY breaking, the potential along a flat direction is given 

by [Hiiniiiaiiii 



U {(j)) = m\ 



|2 



+ ... (30) 



where M is the scale at which ~ ^^3/2 is defined, typically the GUT scale and denote 
non-renormalizable terms with M-p\ suppressed couplings. K encodes the running of the mass 
of the flat direction and is determined by the RG equations (which we discuss at greater length 
in section ISJ, therefore it depends only on the low energy theory below the Planck scale. 

As is well known (see e.g. jjlOj) and rederived in appendix^ the above quantum corrections 
to the flat directions result with either positive or negative pressure, depending on the sign of 
K, 

P K 

If negative, the condensate at the end of inflation fragments and forms Q-balls [5l ll2l05lll3lll4j . 
For the models we are interested in, this process happens relatively quickly after inflation, when 
the VEV of the modulus is of order the Planck scale. Since the charge of Q-balls is proportional 
to the square of the VEV, these large charged Q-balls can be called Q-giants which decay late, 
producing a low reheat temperature (and in some cases, too much baryon number). Thus 
negative pressure flat directions result with unacceptable cosmology and therefore cannot play 
the role of the inflationary moduli. 

After reviewing some known properties of Q-balls in subsection 14.11 we show in 14.21 that 
the condensate fragments into Q-giants (whose statistics we compute) and a non-thermal dis- 
tribution of (f) particles. Section discusses the properties of this gas. 

4.1 Q-Balls Properties 

Let us briefly summarize the basics of non-top ological solitons, otherwise known as Q-balls [Sill] 
(for a review see e.g. /14_ and references therein). Consider a complex scalar field with a U{\) 
invariant potential ?7(|i;^|), 

C = \d^(t^\''-Um. (32) 
The conserved charge is given hj Q = i j d^x{cj)^ (f) — ct>^(t))- We choose an ansatz of the form, 

= 0(r)e^-*, (33) 
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under which the equation of motion (EOM), charge and mass become, 



<f(t> ^ 2# /^-^ (34) 
dr^ r dr ddP' ' 



Q = 2uj j d^x(i?{r), (35) 

Mq = j d^X [{drCP)^ + Uesi^P)] + ioQ, (36) 

with f/ofr(</') = U{(j)) — uj'^cjP'. This EOM describes a particle moving in an inverted effective 
potential with friction oc We seek a solution which minimizes the energy and with boundary 
conditions, 0'(O) = (p'{oo) = 0. A stable solution necessarily exists if 2U{(j)o)/(pQ has a minimum 
away from the origin. This condition holds for the potential (|30|) for negative K, which is the 
case we are dealing with here (postponing positive K to section IH)). 

Here and below, we assume \K\ log |(/)p/M^ <^ 1. This allows us to use the 1-loop approx- 
imation for the potential without having to sum over large logs. We will later see that given K 
of order 10~^ this assumption is easily justified, at least at the onset of Q-ball formation. 

The Q-ball solution relevant in gravity mediated SUSY breaking is the thick wall Q-ball 
solution The ansatz for this solution 



„2 



^{r) = ^oe ^ (37) 

with (f)Q <C Mpi (so that we can ignore non-renormalizable terms). Plugging this ansatz in the 
EOM (jSH) we get. 

Comparing terms one finds 



R ^ y^"v'' ^^^^ 

~ rnl(l-2K + Klog^Y (40) 



and using eqs. p5]) and (|36]) . 



= [^y' 2u>4,lR^ (41) 
Mq - 3g)''%?fl+(2-lA- + 2A-l„g|,)(|)"%„>Sj;3 (42) 
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where we used to ~ and 1/R m^. 

Two remarks are now in order. First, note that eq. ()39() requires that the Hubble radius 
is much larger than m'^^ for Q-balls to form. In particular, H^^/R ^ 1 which leads to 
(j)o/Mpi ^ (3-R'/2)^/^. This requirement is consistent with neglecting the non-renormalizable 
terms. A second remark is that stability is guaranteed by the fact that K < 0. This ensures 
that the effective mass of (p increases as the VEV of (j) decreases and therefore 

^^m,(l-lK- '-Klo, ^) < ^im^ ^ZT^. (43) 

In particular, at 0o = 0, further light degrees of freedom must be reintroduced (recall that the 
above potential is obtained by integrating out heavy particles). Therefore, (pQ in the logs should 
be replaced by the low energy scale which can be taken to be ~ 10^ GeV. Since K < 0, 
this implies that the mass of the free particle [given by the right hand side of Eq. H43|) ] is much 
larger than the mass per charge of the Q-ball, rendering these Q-balls stable. 



4.2 Formation of Q-Giants and Gas 



To predict the size and number of Q-giants, we will make two assumptions which are physically 
reasonable and generic. The first is that, assuming no fine tuning, the orbit of the flat direction 
is homogeneous on a Hubble patch with non- negligible U{1) charge ^H]. The second assumption 
is that in all the relevant length scales, initially, Scpi/cpi ~ H/Mpi ~ 10^"*^^ when (pi is of order 
Mpi and the dynamics we are focusing on begins.^ 

When the condensate fragments, not all particles in the condensate make it into the Q- 
giants. We will refer to the remnant particles as the (p gas. 

Q-Giants: The charge (and therefore size) of these Q-balls is determined by the time at 
which the first perturbation goes non- linear 114j , 

Ho ~ a = - log{5^i^i) 44 

a 

Given that (p dominates the energy density, we can estimate its value when the Q-balls are 
formed, 

, V3HoMpi Vl2\K\Mpi 

(po = ~ (45) 

rUfp a 

^Strictly speaking, the size of the fluctuation might grow somewhat by the time the held moves from a VEV 
> Mpi to the regime of vahdity of this potential. This growth is expected to change the estimate of H at the 
time of Q-ball formation by a multipHcative factor of order 1. In addition, such a process can only serve to 
increase H which will result in an increase of the Q-giants size which will only make our bounds more stringent. 
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and using eqs. (|39jl and (|41j) one obtains, 

^ - |i^|3/2 ^2 - ll03GevJ Vl0-2y V30J ■ ^^''^ 

The validity of this estimate was verified numerically in IT^ . Due to the large initial value 
of (/), these Q-balls store an extraordinary amount of charge, hence the name Q-giants. 

We stress that the (p field that we are using cannot serve as an Affleck-Dine (AD) field 
[371 I17j . The reason is that it dominates the energy density (since it has the largest VEV 
possible in this model) and hence if it were the AD field, it would have generated a very large 
baryon asymmetry. Thus the baryon number must be generated either by another AD field or 
by a variant of EW-baryogenesis [381 1391 1401 1411 1421 1431 144j if the radiation temperature to be 
calculated below is above the electroweak (EW) scale. 

(/>-Gas: The charge in the (p condensate is maximal if the VEV moves in a circular orbit. 
In this case, all the charge makes its way efficiently into Q-giants. In the generic situation 
we are discussing, the orbit is not circular. We show below that in this case, in addition to 
Q-giants, a distribution of low momenta (non-coherent) particles appears. We will refer to 
this as the (p-gas. The existence of such a component is also supported in numerical simulations 
of condensate fragmentation . The critical features of this gas are discussed in section 14.31 

4.2.1 Efficiency of Q-Giants and Gas Formation 

Our initial configuration is that of a condensate, for which the perturbation on some scale grow 
to become non-linear. Once the non-linear regime is formed, Q-giants form very rapidly, as does 
the (p gas. Immediately after fragmentation the universe contains a mixture of the following 
components: 

1. Energy in Q-giants - an order one fraction of the charge initially goes into Q-giants. 

2. A gas of (p particles outside the Q-giants - the remaining energy goes initially into this 
gas. 

3. Light MSSM particles which result from the decay of the (p particles - initially there are 
none of those. 
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The physics of the fragmentation process is the following. The length scale that becomes 
non-linear first is of order 

i.fl-..(^y'^,. (47) 

where a is the expansion factor. Once this mode enters the non-linear regime Q-giants form 
quickly, roughly one per the length scale H47|) . (j) quanta that do not make it into Q-giants 
appear into a non-coherent gas of (p (and anti-0) particles. Since the entire process occurs for 
low momenta - of order up to some powers of if - we expect that the (j) gas will contain 
only low momenta. This picture is also verified numerically in J^l- Furthermore, the same 
simulations suggest that most of the charge finds its way into the Q-giants rather than into the 
(j) gas.^ 

It is also easy to see that the amount of anti-Q-giants is very small. Indeed, recall that the 
condensate is rotating in one direction with a charge close to the maximal allowed. To get a 
Q-giant with the opposite charge, the condensate in a large piece of spacetime has to coherently 
start rotating the other way, which is suppressed. It is also unlikely that a Q-giant will form 
out of a gas of particles after the condensate fragmented. This is from phase space arguments: 
a non-coherent collection of <j) quanta will not be able to assemble itself into a coherent Q-giant 
(or anti-Q-giant). 

The total energy density in the (p condensate at the time of fragmentation is of order 
Ptotai = "T-^^o- '^^^ energy per charge in the condensate is given by ptotai/Q — rn^/{l — e) 
where e denotes the eccentricity of the orbit of (j) and is therefore less than unity. On the other 
hand once Q-giants have formed, they fulfill Mq/Q ~ m<^ [see eqs. H4ip . Thus since 

most of the charge ends up inside the Q-giants, the energy and charge density of the system's 
constituents after the Q-giant formation is 

Qgas/V ~ Pgas ^ e"i^0O 

(48) 

Q/V ~ (1 - e)m^02 PQ^{1- ()ml<pl 

Finally, for consistency of this picture we need to show that the mode-mode interactions 
in the (j) gas are small enough such that the gas does not thermalize but rather remains at low 



®An intuitive argument may be the following. The VEV of the field inside the Q-giant moves in a circular 
orbit in field space. This means that regions in spacetime in which the charge-density/energy-density ratio is 
maximal have an easier time to fall into Q-giants. The gas originate from regions of space in which the charge- 
density /energy- density ratio is low - such regions tend to be less coherent - but this also implies that the gas is 
characterized by a low charge-density /energy-density ratio. 
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momenta until it decays. Indeed, there are only two kinds of interactions which are of impor- 
tance to us: interactions arising from the log term in (|3flj) and non-renormalizable interactions 
coming from corrections to the Kahler potential. For the latter, as discussed in section a 
n-point interaction is suppressed by m^i^/M^^'^ . Thus the rate for momentum transfer can be 
estimated to be 

I K ^^TTV^ 

r;^--!^- - n^M < -^72^, (49) 

where a^f, is the cross-section for these interactions and the \K\'^ factor comes from the suppres- 
sion in the number density of the gas. Here we have estimated Fmm for n = 4. Higher n would 
result with lower rates. 

The log interactions are more subtle. As we discuss at length in section 14.31 while {(j)) 
vanishes outside the Q-giants, the states which interact with remain massive due to non- 
vanishing {(p^ip). Thus the log potential is still valid outside the Q-giants and result with 
interactions suppressed by {(1)^(1)) ■ Thus one finds 

\K\ '^Tfl'^ 

MM ~ ^^t0) ^^^^ 

which is faster than the scattering rate eq. H49() and is of order the decay rate for (p (see 
section In fact, we expect these scattering to be slower than the estimate above since 2-2 
scattering (for which the estimation was done) cannot result with an increase of the average 
momenta while higher n-point scattering is further suppressed. We thus conclude that the gas 
outside the Q-giants is formed out of equilibrium and stays that way until it decays. 

Finally, light MSSM particles are absent until the gas decays. We compute the gas decay 
rate in section In section l5?2] we show that the temperature of the MSSM cannot dissociate 
the Q-giants. 



4.2.2 Initial Velocity and Charge Distribution in Q-Giants 

In section we briefly investigate the consequence of Q-giant collisions. In order to do that, 
we must estimate the distribution of velocity and charge of the Q-giants. We do that here. 



Let us denote a given Q-giant by its center q and radius R. The momentum current is 
then found to be 



P\q,R) 



V{q,R) 



(51) 
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where V{q, R) is the volume of the Q-giant and are co-moving spatial coordinates. To extract 
the distribution of momentum, we follow the evolution of (p until fluctuations become of order 
one, when (p = (pQ. Thus expanding around the homogeneous solution, 

(t){x,t) = (pit) + 5(t){x,t) (52) 

one finds to lowest order in 5(p, 

{SP\q,,R)SP^iqy,R)) ^ / xd^ g- g- {HH^mv)) ■ (53) 

The above can be estimated at the time of Q-giant formation by taking ~ ~ 0o 'ind 
estimating the typical momentum to be of order k ~ R^^. Hence we get 

5P ~ LocplR^ (54) 

which leads to the velocity distribution, 

6P 

6v ~ -— — ~ V l-^l- (55) 

The charge distribution can be found in a similar fashion. Using the expression for Q given 
above, eq. and expanding as in one finds 

{5Q{q,,R)5Q{qy,R)) ^ \cp\^ j d^xd^'ya^ -^-^{5<p\x)5<P{y)) (56) 
from which one can estimate, 

SQ ~ LocplR^ (57) 

and therefore SQ/Q ~ 1. 

Given the above, let us assume an initial Gaussian distribution^'', 

/(*o, Q,v)=N (^^^5^' 

with (Q) given in H41() . Here is the normalization which can be approximated by assuming 
most of the charge is stored in the Q-giants P^ . 

(1 - €)m^4 = jdQQj d\fito,Q,v) (59) 

where e is the initial eccentricity which is of order one (see eq. ()48|1 ) . This gives, 

iV ~ -^i^ ~ (1 - e)m!|i^p/2. (60) 

R'^ ^ 

We will later use this distribution to estimate Q-giant destruction rates. 

^°This assumption is valid in the limit ^/K ^ 1 in which case there are many Q-giants (of order lO/V^) 
within a Hubble radius. 
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4.3 Stability of the Gas 



In the non-relativistic gas of (p particles, the relevant energy per particle are of order ^3/2- 
Hence during the stage of Q-giant formation, one must have, {{dcj))'^ + U{(j))) = which when 
combined with eqs. (|^ and (|15|) gives 

= k^4>l (61) 

where /gas ~ e at the onset of Q-giant formation and is roughly 1 by the time the gas decays 
due to energy loss of Q-giants through collisions (as discussed in section r5.1|) . On the other 
hand (0) ~ between the Q-giants. 

We now argue the following: 

1. All the states which are strongly coupled to <j) heavy - m? ~ (</'^</') - due to this large 
variance in the field and despite the fact that (0) = 0. 

2. While giving large masses to all particles, (j) itself remains light to all order of perturbation 
theory. This is a consequence of super symmetry. 

3. As a result, (j) is stable and decays late, reheating the MSSM particles to temperature 
Tintcr — 'm(j,{MY<\/m(j,Y/^ . This tempertature, as discussed in section cannot dissociate 
the Q-giants, which remain stable and dominate the energy density until they decay much 
later. 

To get a sense for why the above is true, let us work with a specific model. As an example, 
consider the renormalizable interaction^^ 

Wq = X(pX^ + ... (62) 

where A is an order one dimensionless coupling and stands for additional non-renormalizable 
terms. This is the type of interaction we expect in the context of flat directions of the MSSM. 
The scalar potential then takes the form, 

C/scaiar = 4| Aj^ | (/^^ | |X| ^ + |Ap|X|^ + (soft temis) + ... (63) 

and there are Yukawa couplings, 

- ^Yukawa = 2\(j)^ X + XX^^^x (64) 
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(t)t (j) (t)t (j) <j)t (j, 

+ ; I + 

Figure 1: Diagrams contributing to the fermion's two-point functions, 
where are chiral fermions. 

The first term in (|63|) gives a mass to X, m\ = 4|Ap((/)^0), thus making it very heavy. On 
the other hand, one may naively think that ^ x remains hght since it interacts with (p rather 
than 4)^(j). To see that this is not the case, recall that the physical mass for a fermion is not m 
which is charged under chiral rotations, but rather |mp which in that case will be proportional 
to {(t)^4>). Another way to see this is to look in the limit of exact supersymmetry (assuming 
{(l)'^cl)) does not directly break SUSY)-*^^. In that limit X and ^'x must have similar masses and 
therefore the fermion must be massive as well. To be more concrete, let us compute the new 
propagator for "^x in the presence of non- vanishing (0^0). Summing over the diagrams in Fig. 
^one finds the corrected propagator, 

,2 



Indeed, the pole is at 4|A| which is just mj^ . 

One may also worry that thermal and plasma effects induce a large splitting between the 
masses of X and ^'x which contributes to the mass of oc — 'rn\_^. Again, this can 

be argued not to occur. Indeed, in the super symmetric limit, as was argued above, there is 
no splitting between X and ^ x and therefore no contribution the m,^. This can also be seen, 
if one notices that the contribution of the bosonic diagrams to in the presence of (0^0) 
(Fig. cancel the contribution from the fermionic diagrams (Fig. Eb). Any supersymmetric 
breaking effect, arise either from soft breaking terms which are proportional to ~ m3/2 or 
from plasma effects (recall that the gas is not in thermal equilibrium and therefore one cannot 



assign a temperature to it) which must be of order k R ^ = 's/\K\/2m^. Thus any splitting 
must also be of order m3/2 and therefore (p remains light. 

Having established a consistent picture, we may now ask when does the (p gas decay. Since 



^^Of course, if no renormalizable interactions exist and only A/pi suppressed ones, the decay rate of the gas 
is very low 

^■^It is important to distinguish between finite density effects (from and finite temperature SUSY breaking 

effects. 
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Figure 2: Contributions of bosons and fermions to m^. In the supersymmetric limit, the bosonic 
contribution cancels the fermionic contribution exactly. 

it is coupled to light degrees of freedom only through heavy particles of mass ~ {4>^(p), its decay 
rate can be bounded by 

The energy densities in the gas and in the Q-giants are comparable and redshifts like a~^. Thus 
we may assume the gas dominates the energy density and compute the intermediate reheating 
process by solving the equation for the energy density in the radiation, 

Py + AHpy = T^p^. (67) 

Using @ and (|66() and taking H = 2/(3t) for this matter dominated epoch, one can estimate, 

3 /gas 5, /no\ 

This solution is valid until radiation becomes dominant at H = T,^. At this stage, the reheating 
ends and the gas decays completely, producing an intermediate temperature given by 

Tintcr ^ 9~.^"m, (^^^) - 10^ GeV (69) 
with the gas decaying when, 

2/3 



decay V8vr/g4m^y 

Note that unlike the usual inflaton scenario, this intermediate reheating does not start at high 
temperature and therefore there is no danger to the stability of the Q-giants at this stage. 
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From that point onwards, radiation surrounds the Q-giants. Initially it dominates the en- 
ergy density but very quickly (once the expansion factor grows by two orders of magnitude) it 
becomes sub-dominant with Q-giants dominating until they decay. Note that taking into ac- 
count sphalerons in the case where the charge in the giants corresponds to B-|-L, is not expected 
to change the final B-l-L charge since the gas outside the Q-giants stores very little charge. Be- 
low we further investigate the dynamics of the Q-giants and show that the temperature of 
radiation does not dissociate the Q-giants. 

5 Evolution of Q- Giants and Reheating 

The evolution and stability of Q-balls has been widely investigated j451 |H1 111! 1141 146j . There 
are many effects which can influence the time of decay and the reheat temperature. We will 
consider the most relevant effects: fragmentation of Q-giants through collisions, dissociation 
due to external thermal plasma, and evaporation to bosons and fermions. 

In addition, there is potentially a new effect which is the decay of Q-giants when immersed 
in the non-thermal gas of low momenta (p particle. However, as we saw above, in cases where 
the Q-giants are large the gas outside has large ((/<V) ^ fgas4>o which suppresses interactions 
between the Q-giants and the <p gas at the surface (as well as with light particles). 

5.1 Q-Giant Collisions 

Q-giants may collide between themselves leading to charge and energy transfer and may be the 
cause of Q-giants disintegration. Collisions are especially important since the reheat temper- 
ature is inversely proportional to a power of the charge (see equation H92|)). and hence if the 
Q-giants fragment to Q-balls with lower charge, the reheat temperature will be higher. We now 
show that Q-giant collisions only mildly reduce the number of Q-giants with large charge. 

The only long range Q-giant /Q-giant interaction is gravitational. In addition they interact 
through direct collisions and interact with the gas around them. Here we focus on the effects 
of collisions. We expect that the cross-section for collision is of order the size of the Q-giant, 
which is verified in numerical simulations ^3B^. More precisely, in general collisions are divided 
into three types: fusion, charge exchange and elastic. The cross-section for each is different 
with charge exchange being the largest, of order four times the geometric cross-section^^. To 
^■^This is the largest cross-section for a process that can reduce the charge. A larger cross-section can occur 
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see how collisions reduce the number of Q-giants, let us define 

/■3(Q>/2 r 

ngit)^ / dQ' / d\f{t,Q\v) (71) 

J{Q)/2 J 

with /(t, Q, v) is the distribution function in velocity v and charge Q of the Q-giants (as defined 
in ((2H1))- We assume that whenever there is a collision between Q-balls, they fragment com- 
pletely. With this assumption we are underestimating the number of Q-balls, and in particular 
the large Q-giants (hence obtaining an upper bound on the reheat temperature). 

Within this approximation, the equation of motion for nq is given by, 

dnq 



dt 



+ ZHnq = -n^qav (72) 



with a estimated to be of order the geometric cross-section a ~ vri?^, which is approximately 
constant. In the term on the right hand side we have neglected some order one coefficient which 
is related to the fact that the decay width is proportional to the total number of Q-balls and 
not just the Q-giants defined in H71() . As was shown in section [4.2.21 we may take the initial 
velocity, which redshifts like a~^, to be of order vq ~ y^|/C|. Thus using 1)441) . one finds 



+ J J =nq[-) e^Ma) (73) 

where Hq ~ €m^\K\^^'^ is the initial number density for Q-giants with charge of order (Q). 
We see that collisions become less effective as the universe expands. This is easily understood 
since, as the number density of the Q-giants drops down, collisions become more rare. Taking 
again a ~ 30, we learn that collisions reduce the number of Q-giants by roughly two orders of 
magnitude, namely 

ecou(a^oo) > 10-2. (74) 

Thus even if we assume that collisions break Q-giants only into (p gas (without small Q- 
balls), we see that still a few percent of the energy density is left in the Q-giants before the (j) 
gas decays. Thus we may take, 

/gas 1 - ecoii(a oo), (75) 

[see eq. H61() ]. 



for fusion, but this will make the Q-giants larger. The cross-section for fusion strongly depends on the phase 
correlation between colliding Q-giants. Q-giants with a similar phase will tend to merge, while the fusion cross- 
section for randomly distributed phases is highly suppressed. 
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Since we have assumed that the coUisions cause the Q-giants to fragment, this is a lower 
bound on the number of Q-giants. Two additional mechanisms change the survival rate of the 
Q-giants upwards by a small amount: 

1. The field theory binding energy of the Q-balls is of order |Er|MQ, where Mq is the mass 
of the Q-ball. The kinetic energy is initially of order Mqv^ /2 ~ MqK. The kinetic energy, 
however, redshifts rapidly and very quickly there is not enough kinetic energy to overcome the 
binding energy. This will increase the survival rate of the Q-giants by a factor of order 1. 

2. In the simulations [ISl Q-balls often do not fragment upon collision. In fact, for the 
Q-giants the gravitational binding energy makes them even more resistant to fragmentations 
since the gravitational energy, 



mlQ^ ^ 104 \K^Mf 



is proportional to Q'^ and is similar in order to the binding energy of the field in flat space (see 
(jlHl)) which is 

^binding ^ K m^Q ~ ^ ^. 77 

Egrav is proportional to Q^, which means that large Q-giants are energetically favored over 
several smaller Q-giants of the same charge. However, for a ~ 30, this effect is rather small. 

The conclusion is that by the time the collisions end - shortly after the Q-giants form - 
about 1% of the total charge and energy is in non-relativistic Q-giants, and the remaining charge 
and energy is in the (f) gas. Hence, when the gas decays, there is a short period of radiation 
domination, before the Q-giants become dominant again. As we show in the following sections, 
given the intermediate temperature, eq. H69jl . there is enough time between the decay of the 
gas and the decay of the Q-giants for the Q-giants to take over and dominate of the energy 
density. 



5.2 Dissociation of Q-Giants 



After the 0-gas decayed, the universe can be thought of as two subsystems composed of Q- 
giants and radiation. Because the gas never reached thermal equilibrium before decaying, the 
two systems are not in thermal equilibrium. Nevertheless, the effects of thermal radiation may 
be important and must be taken into account. We now compute a bound on the temperature 
outside the Q-giants below which, the Q-giants dissociate. To do so, we briefly repeat the 
considerations presented in [11) . 
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Thermal particles can penetrate Q-giants thereby transferring energy and possibly causing 
dissociation. The condition to evade dissociation is that the excess energy per unit charge 
(which was not radiated away) is smaller than the binding energy per unit charge 6m ~ \K\m(p: 

< 6m (78) 

5tr 



Q 

where 6tr is the time scale over which the Q-giants radiates its excess energy away. The flux of 
particles hitting the Q-giant is 

r = ^A^rrlT' (79) 

where Vst is the stopping radius of the Q-giant defined by g4){rst) — 3T with g ~ 0(1) being 
the coupling constant of (j) to radiation and g^, is the number of relativistic degrees of freedom 
coupled to the Q-giant. Taking the energy per thermal particle transferred to the Q-giant be 
7Tr, we get 

dMQ ^ 47T<7*(T)o^4 



dt vr 
Tst can be easily calculated by using the ansatz (|37|) . 



Thus, using (|39|) for R we get, 



MqT = ' ritT\ (80) 

(81) 



dMQ 8^Tg*iT)l3^ 
at TT\K\m, 

which, from the inequality (f78|) and using 6tr ~ l/m^ one obtains 

where we have taken 7t ~ 1 and /? ~ 4. 

In lU this inequality is used as a bound on the reheat temperature. Here, we use it as 
a bound on the maximal temperature outside the Q-giants. It is therefore immediately clear 
that the radiation after the decay of the gas cannot dissociate the Q-giants. Below, we also 
investigate the reheating process during the decay of the Q-giants themselves. We show that 
dissociation does not play a role at that stage either. 

5.3 Reheating 

In general Q-balls can decay into fermions and bosons through the surface where these particles 
are light |45l lllj . The decay into fermions was found to be bounded by j45j 



dtdA - 1927r2 ' 
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with = u;-^ — 771^(1 + K) ~ 3m^|iC|, whereas the decay into bosons, as was shown in is 
enhanced by a factor fs — 170/|i^|^/^.^^ From the above, one obtains an upper bound on the 
decay width of the Q-giants^^ 



Let us now evaluate the temperature outside during and at the end of the reheating process. 
To do so, we solve the Boltzmann equations, 

-^+^Hp^ = / dg/(i,Q)rg-M(Q), (86) 

^I^ + 2,Hf{t,Q) = QT^-^^IMl. (87) 

where /(t, Q) = f d'^vf{t, Q, v) is the time-dependent density distribution function for a Q-ball 
with charge Q, p-y is the energy density stored in radiation and M{Q) ~ Qm^ is the mass of 
the Q-ball (see eq. (02) )• The last term in eq. ()87|) is a transport term in Q-space, which takes 
into account the fact that as the Q-giant decays, it loses its charge gradually. The solution for 
(|87|) is given by, 

/(t,Q) = '/(to,Q + fg"^(t-to) 

where T'q'^^ = QVq"^^ ~ is independent of Q. Thus, f{t,Q) is simply the initial Gaussian 
distribution with a shift of its central value, (Q) (Q) — rQ^^(t — Iq). This solution is easy 
to understand. For small t, t — Iq <^ {Q) /^q'^^ = (^'{Qf)~^^ distribution is approximately 
constant and the Q-giants redshift like matter. At H^^ = t — to = (r^^)~^ the center of the 
distribution becomes negative where the solution is no longer valid. At this point, the Q-giants 
decay instantaneously transferring the energy to radiation. 

Next, we solve for the energy density in radiation. The solution to H86() is 

= (—) ' fdt' (^)' r'^^\Qf{t\Q)T^^-M{Q). (89) 

\ '^O / ^ ti V '^O / J Qmin 



Using the above approximation of constant /(t, Q), we get for t — to <^ (T'^q^) ^, 



3(1 - e)m^|i^|3/2 
Pi- — 



5Ho \ao 



a 



5/2 

- 1 -1 



a \ ' 



(90) 



^^Note that the actual decay rate for the Q-giant scenarios is much smaller initially until the gas between the 
Q-giants decay. 

^^The reader should note that finite temperature decay, as was computed in I47| is not relevant in our scenario 
since the Q-giants are very weakly coupled to the light particles. 
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where we took Qmin = and Qmax = c>o, and Hq is given in (|44|) . Thus at the transient time, 
radiation reheats to a maximal temperature, 

and redshifts hke a"'^/^ until t ~ ^F^.^^^ when the Q-giants decay completely, reheating the 
universe to 

Once the Q-giants decayed, radiation redshifts like as usual. 

As was shown in 14.31 the intermediate reheat temperature from the decay of the gas is some 
nine orders of magnitude above the reheat temperature of the Q-giants. This, as we advocated 
above, is sufficient to allow the Q-giants to regain domination over the energy density, lost due to 
collisions. Furthermore, given the bound on dissociation, eq. ()83() . we thus find that during the 
reheating process, Q-giants are not dissociated and the final reheating temperature is very low 
due to their large charge. Hence, the Q-giants dominate the energy density too long, leading to 
a too low reheating temperature for a successful conventional big bang nucleosynthesis required 
in the standard cosmological scenario. We are therefore led to conclude that negative pressure 
flat directions cannot play a role in modular inflation models. 



6 Positive Pressure Flat Directions and Primordial Black Holes 



We now turn to flat directions which, due to RG running, generate positive pressure, namely, 
K > 0. Such flat directions do not fragment to form Q-balls, but rather generate a finite 
Jeans scale which prevents the formation of primordial black hole after inflation ^S]. Indeed 
since typically the modulus that drives inflation is weakly coupled, there is a long stage of 
matter domination after inflation which can generate large density perturbations and lead to 
the creation of BHs ^Hj. It is interesting that within the MSSM, the only direction which is 
not charged under B-L and which has positive pressure is the HuHd direction. 

Consider again the potential (|3Ujl . only this time with positive K. As is shown in Appendix 
El the equation of state is approximately given by, 

K K 
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and the Jeans scale is therefore 

Xj = ^H-\ (94) 

At scales smaller than Aj pressure prevents growth of perturbations and such perturbations 
oscillate. Thus whether or not BHs can form, depends on the value of K at the time density 
perturbations becomes non-linear. K is given by 

1 dml 



K 



m 



2 



dt 



(95) 

i=log{|</.p/M2) 



with t = log(/iV^^) and can therefore be extracted from the RGE's of a given model. In what 
follows, we assume the MSSM up to the GUT scale. Our strategy is to find the Jean's scale at 
the time where 6p/ p is of order one and compare it with the scale of the fluctuations. If the 
Jean's scale is larger we are assured that no Q-balls or BHs can form. 

Assuming 5p/ p ~ 5 x 10^^ for scales of order the horizon at the end of inflation, and using 
the fact that during matter domination, perturbations grow like a while redshift like a~^/^, 
we find 



^^ = I^I.P/P=i = (5 , io")-3/2 - lO^'GeV. (96) 

On the other hand, a density fluctuation of size Ainf ^ at the end of inflation becomes 
when 5p/ p = 1, 

A|5p/p=i = a-^'^H-^ ~ 2.2 X lO^^i?^^ (97) 

Thus to a-priori eliminate the possibility of BH formation we must look for a flat direction of 
the MSSM with K/2 Z 2.2 x lO^^. 

The RGE's for the MSSM are given in jl^. For the particular case of HuH^ flat direction, 
one finds^^ K/2 ~ 2.3 x 10^^. We can therefore conclude that no primordial BHs can form 
if the HuHd direction plays the role of the modulus in inflation. Moreover, it is the only 
possible direction in the MSSM. Before structure can form, the modulus must decay and reheat 
the MSSM. After that, structure formation can follow as usual, in the radiation dominated 
universe. 

Note that we have imposed a very strong restriction, which is fortunately fulfilled in the 
MSSM. It could well be that 6p/p ~ 1 is permissible, in which the universe undergoes a period 
of structure formation by (j) particles, similar to cold dark gas. It will be interesting to explore 
this possibility further. 

^®We thank Guy Englehard for enabling us to use his Mathematica code for the RGE's. 
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7 Discussion and Caveats 



In this section, we briefly detail what would go wrong with cosmological scenarios in which 
the negative pressure flat directions survive and argue that thermal inflation [231 1241 125j 
may change the conclusions. However, we present reasons for the belief that invoking thermal 
inflation in the present context requires fine tuning. 

With Q-giants surviving a long time, big bang nucleosynthesis would be severely disrupted 
since the historical sequence of radiation domination (RD) and matter domination (MD) periods 
will be modified from the standard scenario of 

inflation + MD ^ RD(BBN) ^ MD ^ today (98) 

to 

inflation + MD ^ RD ^ MD(BBN) ^ RD ^ MD ^ today (99) 

or 

inflation + MD ^ RD ^ thermal inflation + MD ^ RD (BBN) ^ MD ^ RD ^ MD ^ today. 

(100) 

Eq. (|99() depicts the situation in which there is no thermal inflation, while Eq. HlOOj) depicts the 
situation in which there is thermal inflation to dilute possibly too large baryon number (note 
that we have included any possible additional period of matter domination between thermal 
inflation and radiation domination following from flaton decay). The label "BBN" marks the 
period during which nucleosynthesis occurs. A successful nucleosynthesis is impossible by the 
time of last RD periods labeled in Eqs. (|99|) and HlOOj) since the maximal reheating temperature 
of 0.1 MeV during those time periods is smaller than ~ 10 MeV required for chemical equilibrium 
initial conditions or ~ 1 MeV when neutrons drop out of equilibrium [SUI . 

Clearly, without thermal inflation BBN occurs during matter domination and thus the 
existence of Q-giants is excluded. One can also show that the baryon-to-entropy ratio during 
BBN in that case does not allow for a successful production of the light elements. We hence 
consider the scenario depicted by Eq. (|100() in which one invokes a short period of thermal 
inflation [231 1241 l25j . Assume that we have the flaton, (/? (not to be confused with 0), at some 
VEV, M, initially at its true vacuum. Once the gas decays, reheating the outer region of 
the Q-giants to a maximal temperature of Tbegin = Winter = 10^ GeV, thermal corrections 
drive i-p to the origin, and thermal inflation begins at roughly that temperature. It ends when 
^'^We thank Ewan Stewart for drawing our attention to this issue. 
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Tend = n^ip — 10^ GeV and thus there is at most 5 e-folds of thermal inflation, which can reduce 
the energy density in the Q-giants by 10~^. 

The energy density during thermal inflation is of order Vq = m^M^. Since the flaton stores 
no energy when it sits in the true minimum, the total energy density once the gas decays is of 
order g*T]nter therefore Vq must be smaller, giving. 

The above can be understood in a different way, by noting that at M, ip is coupled to radiation 
through massive particles with mass M. Thus corrections to its mass are expected to be of 
order T'^/M'^ and therefore, for these corrections to drive it to zero, H1U1|1 must hold. 

Within the context of gravity mediated supersymmetry breaking one can argue that such 
a requirement from thermal inflation is difficult to achieve ^5]. Indeed, the flaton is expected 
to have soixLG 'Zjji symniGtry {ji ^ 2) with. a. supGrpotGnticil, 

W = ^-^ + ... 102) 

Taking into account soft breaking mass —rn^i^ip^ and A-terms 1^3/2^ one finds the minimum 
of the potential to be at 

M = {ip)^ Mpi ^-j^j . (103) 

Hence the only value of M which is compatible with 1)101(1 is M = 711,3/2 — 10^ GeV for n = 3. 
However, in such a case, thermal inflation begins only at Tbegin = 10^ = Tend and so there is no 
time for thermal inflation to dilute the Q-giants. Thus invoking a period of thermal inflation 
in this scenario requires fine tuning. Nevertheless, for completeness, let us suppose that one 
is able to produce such an epoch and see what are the requirements for thermal inflation to 
sufficiently dilute the Q-giants. 

At the end of thermal inflation, assuming perturbative decay of the flaton, the energy 
density in radiation is of order p^^^~'^'^'^'^^ {Tfi^.^) = ^m^M'^ where 7 is a factor coming from the 
volume dilution during the coherent oscillation of the flaton and 10^^ GeV < Tgat < 10'^ GeV is 
the flaton's reheating temperature (note that the lower bound on the flaton decay temperature 
comes from the requirement of successful BBN). On the other hand, the energy density in 
Q-giants is (including collisions - see section [^?T|) . 

PQ(Tfiat) ^ 10-^5* (Tbegin)r4gin- (104) 
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After the flaton decays, in principle, there could occur a baryogenesis scenario and then BBN can 
commence afterwards. Now, for BBN not to be disrupted, we must have radiation domination 
during that time. This imphes, 

The lower bound decreases for smaller Tfjat because the flaton decay radiation energy density 
dilutes faster with the scale factor than when the flaton energy density is stored in the form 
of coherent oscillations. In the limit that the flaton never decays but remains in coherent 
oscillations, since the Q-giants energy density dilutes the same way as that of the flaton, they 
never come to dominate the energy density. Eqs. HlOlj) and H105|) mean that in principle, 
thermal inflation can allow those flat directions not carrying B to survive. 

If the Q-giants carry B charge^^ the situation is more constrained. In order not to produce 
too much baryon number, the Q-giants must never dominate the energy density. In such a case 
it is simple to compute the baryon asymmetry at the time of Q-giant decay, Treheat — 10~^. 
Indeed, at reheat, 

PqiTreheat) — PQ(^flat 

while, 

P7(^flat) ^ im^M^ 

Hence using ub/s ~ {pQ/'m<j,)/p^^'^ and by demanding that the baryon-to-entropy ratio so 
produced is smaller than the measured 10"^'', one finds 

M > f r^/y ^flat y^V^y/^G,v. (108) 

V103GeV/ VlO^GeVy V200/ ^ ^ 

Thus thermal inflation can, in principle, dilute the baryon number produced by Q-giants decay, 
if its reheat temperature is between 10~^ GeV < Tfjat < 10 GeV. 

Finally, consider the flat direction with positive pressure such that no Q-giants form. In 
that case, thermal inflation can acceptably dilute the baryon number density if 

M > IQio ( -^^) ( ( 3*siTnn)\ ^ ^ 

V 105 GeV y Vio^Gevy V 200 y ^ ^ 

^*Note that phenomenologically, it is more appropriate to classify the scenarios according to baryon number 
since the observational constraints on lepton number is very weak |51| . 



V ^flat / 



^reheat \ 
^fiat / 



(107) 
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where Trh < lO^GeV is the reheating temperature after inflation ends^^ and just as before, Tflat 
corresponds to the temperature at which the flaton decays. For sufficiently small Trh and Tgat 
one can dilute the produced baryon number to an acceptable value. Hence, the flat direction 
HuL for which K > can survive with sufficient engineering of thermal inflation. (Note that if 
the reheating temperature is sufficiently low, sphalerons may be sufficiently suppressed to forbid 
the conversion of the Icpton number into baryon number. In that case, no further dilution from 
thermal inflation may be necessary.) 

To conclude this section, we point out that while thermal inflation can, for a narrow win- 
dow of the reheating temperature, dilute the Q-giants sufficiently thus evading the constraints 
considered in the previous sections, such a theory requires fine tuning and so seems to be im- 
probable. The same conclusion holds for positive pressure flat directions. Finally, even for the 
flat directions that are generically "excluded," extremely unlikely fine tuning of initial condi- 
tions (such that no charge is produced) may allow them to be viable. 

8 Conclusions 

It has long been realized that negative pressure flat directions fragment into Q-balls which are 
rather stable and may survive the EW phase transition. In this paper, wc have considered the 
evolution of such flat directions that arc assumed to play a role in inflation, and which therefore 
have VEVs of order the Planck scale at the end of inflation and which dominate the energy 
density at that time. Such flat directions are found, for example, in hybrid inflation models. We 
have shown that within the framework of gravity mediated supersymmetry breaking, these flat 
directions fragment into very large "Q-giants" which survive nucleosynthesis and are therefore 
excluded cosmologically. 

The small decay width of the Q-giants is directly related to its large charge. However, 

thermal effects may break and dissociate the Q-giants long before they decay. Indeed, a new 

effect which was not considered in previous works is the formation and decay of a gas outside 

the Q-giants. During the fragmentation of a flat direction, 0, not only Q-giants are formed, 

but also (p particles in the form of a non-thermal gas. Although, the majority of the charge 

is stored inside the Q-giants, there is an order one fraction of energy density stored in the 

gas. While outside the Q-giants {(p) vanishes, the naive expectation, that the ^-gas decays very 

^^Note that the finite density effect applies independently of whether there is negative pressure or not. Note 
further that the upper bound of 10® GeV on the temperature of the universe is quite robust in such models where 
the modulus controls the energy density. 
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quickly and reheats the universe to high temperature which would dissociate the Q-giants, is 
wrong. The reason is that {(t)^ (p) generates large masses to the particles to which it couples, and 
therefore its decay is mediated through these particles with masses of order (0^^). One finds 
the intermediate temperature outside the Q-giants, once the gas decays to be, 

Tinter - lO' GeV, (110) 

which is much less than the dissociation temperature, Tjiss — 10® GeV. Thus the temperature 
outside the Q-giants is never high enough to dissociate them. 

Taking collisions into account, we have also shown that by the time the gas decays it 
dominates the energy density and therefore there is a short period of radiation domination, 
before the Q-giants take control over the energy density until they decay and reheat the universe 
to the temperature 

Treheat < 10"^ GeV (111) 

which is generically too low of a reheating temperature to be compatible with cosmology. 

Hence, without further entropy-releasing mechanisms, negative pressure flat direction sce- 
narios are not viable cosmological scenarios. As an example of entropy-releasing mechanism, 
we have considered thermal inflation and found that such a scenario requires fine tuning. Fur- 
thermore, such period of inflation must be properly engineered, with only a narrow window for 
the reheating temperature that allows these fiat directions to become viable again. 

The only fiat directions that do not suffer from too low reheating temperature without 
thermal inflation are positive pressure flat directions. Hence, within the MSSM, the only 
possible directions are HuH^ and HuL. Since the latter is charged under B-L and if the reheating 
temperature is large enough to produce baryon number, it requires a dilution mechanism like 
thermal inflation (which again requires fine tuning) . Finally, we point out that these directions 
prevent the formation of primordial BHs due to a generation of a finite Jeans scale. 
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A The Virial Theorem and Equation of State 

We present here a simple derivation to the pressure generated by a scalar condensate oscillating 
in the potential H30|) which can also be approximated by 

U = mV^J^- (112) 

This result was first derived in jlUj . 

To do so we consider an alternative version of the virial theorem. Consider the quantity, 

(F) = (|(<^<^t + <^^t)) (113) 

where the average is taken with respect to time. Hence, as long as the motions are bounded, 
we have 



1 ■ 

(F) = hm - / Fdt 



Thus we can write, 





lim l[F(r)-F(0)]=0. (114) 

1 — >oo 1 



(2|</>|+</><^t + 00t) 

(2|0| - 2m2 |c/>|2(l + K + K\n ^)) (115) 



where we have used the potential (|3flj) and ignored the expansion of the universe. This is 
justified assuming the oscillation period is much smaller than . We therefore find, 

= m2(|<^|2(l + K + KlnM!)). (116) 



Next, we have. 



p = (|<^|2 + t/)=m2(|</>|2(2 + K + 2KlnM!)), (117) 
P = {\4^\^-U)=mlK{\<t^\^), (118) 
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which gives us 

P K , , 

- = ur2 ■ 119 

P 2 + K + {2K\cf>\Hn%)/{m 



Similarly, using the potential (|112|1 on finds, 

P K 



(120) 



p 2 + K 
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